Non-coding RNAs in the nervous system

Increasing evidence suggests that the development and function of the nervous system is heavily dependent on RNA editing and the intricate spatiotemporal expression of a wide repertoire of non-coding RNAs, including micro RNAs, small nucleolar RNAs and longer non-coding RNAs. Non-coding RNAs may provide the key to understanding the multi-tiered links between neural development, nervous system function, and neurological diseases.

A long non-coding RNA targets microRNA miR-34a to regulate colon cancer stem cell asymmetric division.

The roles of long non-coding RNAs (lncRNAs) in regulating cancer and stem cells are being increasingly appreciated. Its diverse mechanisms provide the regulatory network with a bigger repertoire to increase complexity. Here we report a novel LncRNA, Lnc34a, that is enriched in colon cancer stem cells (CCSCs) and initiates asymmetric division by directly targeting the microRNA miR-34a to cause its spatial imbalance. Lnc34a recruits Dnmt3a via PHB2 and HDAC1 to methylate and deacetylate the miR-34a promoter simultaneously, hence epigenetically silencing miR-34a expression independent of its upstream regulator, p53. Lnc34a levels affect CCSC self-renewal and colorectal cancer (CRC) growth in xenograft models. Lnc34a is upregulated in late-stage CRCs, contributing to epigenetic miR-34a silencing and CRC proliferation. The fact that lncRNA targets microRNA highlights the regulatory complexity of non-coding RNAs (ncRNAs), which occupy the bulk of the genome.

Dissecting the Functional Impacts of Non-Coding Genetic Variation

A large proportion of the variation in traits between individuals can be attributed to variation in the nucleotide sequence of the genome. The most commonly studied traits in human genetics are related to disease and disease susceptibility. Although scientists have identified genetic causes for over 4,000 monogenic diseases, the underlying mechanisms of many highly prevalent multifactorial inheritance disorders such as diabetes, obesity, and cardiovascular disease remain largely unknown. Identifying genetic mechanisms for complex traits has been challenging because most of the variants are located outside of protein-coding regions, and determining the effects of such non-coding variants remains difficult. In this dissertation, I evaluate the hypothesis that such non-coding variants contribute to human traits and diseases by altering the regulation of genes rather than the sequence of those genes. I will specifically focus on studies to determine the functional impacts of genetic variation associated with two related complex traits: gestational hyperglycemia and fetal adiposity. At the genomic locus associated with maternal hyperglycemia, we found that genetic variation in regulatory elements altered the expression of the HKDC1 gene. Furthermore, we demonstrated that HKDC1 phosphorylates glucose in vitro and in vivo, thus demonstrating that HKDC1 is a fifth human hexokinase gene. At the fetal-adiposity associated locus, we identified variants that likely alter VEPH1 expression in preadipocytes during differentiation. To make such studies of regulatory variation high-throughput and routine, we developed POP-STARR, a novel high throughput reporter assay that can empirically measure the effects of regulatory variants directly from patient DNA. By combining targeted genome capture technologies with STARR-seq, we assayed thousands of haplotypes from 760 individuals in a single experiment. We subsequently used POP-STARR to identify three key features of regulatory variants: that regulatory variants typically have weak effects on gene expression; that the effects of regulatory variants are often coordinated with respect to disease-risk, suggesting a general mechanism by which the weak effects can together have phenotypic impact; and that nucleotide transversions have larger impacts on enhancer activity than transitions. Together, the findings presented here demonstrate successful strategies for determining the regulatory mechanisms underlying genetic associations with human traits and diseases, and value of doing so for driving novel biological discovery.

Investigation of the role of long non-coding RNAs in oncogene induced cellular senescence

Cellular senescence is a stable arrest of cell proliferation induced by several factors such as activated oncogenes, oxidative stress and shortening of telomeres. Senescence acts as a tumour suppression mechanism to halt the progression of cancer. However, senescence may also impact negatively upon tissue regeneration, thus contributing to the effects of ageing. The eukaryotic genome is controlled by various modes of transcriptional and translational regulation. Focus has therefore centred on the role of long non- coding RNAs (lncRNAs) in regulating the genome. Accordingly, understanding how lncRNAs function to regulate the senescent genome is integral to improving our knowledge and understanding of tumour suppression and ageing. Within this study, I set out to investigate the expression of lncRNAs’ expression within models of senescence. Through a custom expression array, I have shown that expression of multiple different lncRNAs is up-regulated and down regulated in IMR90 replicative senescent fibroblasts and oncogene-induced senescent melanocytes. LncRNA expression was determined to be specific to stable senescence-associated cell arrest and predominantly within the nucleus of senescent cells. In order to examine the function of lncRNA expression in senescence, I selected lncRNA transcript ENST0000430998 (lncRNA_98) to focus my investigations upon. LncRNA_98 was robustly upregulated within multiple models of senescence and efficiently depleted using anti-sense oligonucleotide technology. Characterisation and unbiased RNA-sequencing of lncRNA_98 deficient senescent cells highlighted a list of genes that are regulated by lncRNA_98 expression in senescent cells and may regulate aspects of the senescence program. Specifically, the formation of SAHF was impeded upon depletion of lncRNA_98 expression and levels of total pRB protein expression severely decreased. Validation and recapitulation of consequences of pRB depletion was confirmed through lncRNA_98 knock-out cells generated using CRISPR technology. Surprisingly, inhibition of ATM kinase functions permitted the restoration of pRB protein levels within lncRNA_98 deficient cells. I propose that lncRNA_98 antagonizes the ability of ATM kinase to downregulate pRB expression at a post-transcriptional level, thereby potentiating senescence. Furthermore, lncRNA expression was detected within fibroblasts of old individuals and visualised within senescent melanocytes in human benign nevi, a barrier to melanoma progression. Conversely, mining of 337 TCGA primary melanoma data sets highlighted that the lncRNA_98 gene and its expression was lost from a significant proportion of melanoma samples, consistent with lncRNA_98 having a tumour suppressor functions. The data presented in this study illustrates that lncRNA_98 expression has a regulatory role over pRB expression in senescence and may regulate aspects of tumourigenesis and ageing.

FROM GENES TO GENOMES: LOCAL ADAPTATION AND ADAPTIVE POTENTIAL IN TWO ARCTIC SEABIRDS

By investigating the mechanisms underlying the evolution and the maintenance of local adaptations we can help predict how species will adapt to future environmental change. In this thesis I investigate local adaptation and adaptive potential in thick-billed and common murres (Uria lomvia and U. aalge), two arctic seabirds of international conservation concern. Thanks to the recent development of new genomic methods, I address three major themes that are relevant for both the development of evolutionary theory and conservation: 1) the role of gene flow in the origin and maintenance of adaptation; 2) levels and distribution of standing genetic variation, and their contribution to adaptive potential; and 3) the genomic mechanisms maintaining an adaptive dimorphism within a single interbreeding population. First, I review the literature on genomics of local adaptation with gene flow and find that adaptation can be maintained despite gene flow, that gene flow itself can promote adaptation, and that genetic architecture is important in the origin and maintenance of local adaptations. Second, I genotype genome-wide markers and toll-like receptor genes (TLRs) to investigate local adaptation and adaptive potential in thick-billed murres. Thick-billed murres do not show signatures of local adaptation to their breeding grounds, but outlier loci group birds according to their non-breeding distributions, suggesting that selection and/or demographic connectivity in the winter may explain patterns of differentiation in this species. Genetic variation at TLRs does not decrease with increasing latitude as predicted, but tests of selection and measures of genetic diversity suggest differences in local selective regimes at most genes. Thick-billed murres show high levels of standing genetic variation and their adaptive potential will mostly depend on rate and magnitude of environmental change. Finally, I improve and annotate the assembly of the highly heterozygous genome of the thick-billed murre. Using this assembly as a reference, I perform whole genome analyses to investigate the genomic basis of an adaptive dimorphism in Atlantic common murres. I show for the first time that a 60 kb complex copy number variant in a non-coding region maintains differences in plumage and cold adaptation despite high gene flow.

The complete mitochondrial genomes of two common shrimps (Litopenaeus vannamei and Fenneropenaeus chinensis) and their phylogenomic considerations

Complete mitochondrial genomes have proven extremely valuable in helping to understand the evolutionary relationships among metazoans. However, uneven taxon sampling may lead to unclear or even erroneous phylogenetic topologies. The decapod crustaceans are relatively well-sampled, but sampling is still uneven within this group. We have sequenced the mitochondrial genomes of two shrimps Litopenaeus vannamei and Fenneropenaeus chinensis. As seen in other metazoans, the genomes contain a standard set of 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and an AT-rich non-coding region. The gene arrangements are consistent with the pancrustacean ground pattern. Both the pattern of gene rearrangements and phylogenomic analyses using concatenated nucleic acid and amino acid sequences of the 13 mitochondrial protein-coding genes strengthened the support that Caridea and Palinura are primitive members of Pleocyemata. These sequences, in combination with two previously published penaeid mitochondrial genomes, suggest that genera within the family Penaeidae have the following relationship: (((Penaeits + Fenneropenaett.) + Litopeiiaelts) + Marsupenaeus). The analyses of nucleic acid and amino acid sequences of the mitochondrial genomes also strongly support the monophyly of Penaeidae, Brachyura and Pleocyemata. In addition, the analyses of the average Ka/Ks in the 13 mitochondrial protein-coding genes of penaeid shrimps indicated a strong purifying selection within this group.

Measles virus minigenomes encoding two autofluorescent proteins reveal cell-to-cell variation in reporter expression dependent on viral sequences between the transcription units.

Transcription from morbillivirus genomes commences at a single promoter in the 3' non-coding terminus, with the six genes being transcribed sequentially. The 3' and 5' untranslated regions (UTRs) of the genes (mRNA sense), together with the intergenic trinucleotide spacer, comprise the non-coding sequences (NCS) of the virus and contain the conserved gene end and gene start signals, respectively. Bicistronic minigenomes containing transcription units (TUs) encoding autofluorescent reporter proteins separated by measles virus (MV) NCS were used to give a direct estimation of gene expression in single, living cells by assessing the relative amounts of each fluorescent protein in each cell. Initially, five minigenomes containing each of the MV NCS were generated. Assays were developed to determine the amount of each fluorescent protein in cells at both cell population and single-cell levels. This revealed significant variations in gene expression between cells expressing the same NCS-containing minigenome. The minigenome containing the M/F NCS produced significantly lower amounts of fluorescent protein from the second TU (TU2), compared with the other minigenomes. A minigenome with a truncated F 5' UTR had increased expression from TU2. This UTR is 524 nt longer than the other MV 5' UTRs. Insertions into the 5' UTR of the enhanced green fluorescent protein gene in the minigenome containing the N/P NCS showed that specific sequences, rather than just the additional length of F 5' UTR, govern this decreased expression from TU2.

La régulation du gène CYP19A1 dans les cellules de granulosa bovine in vitro

L’oestradiol joue un rôle important dans la reproduction en général, particulièrement dans la croissance folliculaire chez la vache. La production de l’œstradiol nécessite l’expression du gène CYP19A1 suite à la stimulation des cellules de granulosa par l’hormone folliculostimulante (FSH) ou le facteur de croissance insulinique de type 1 (IGF-1). Chez la vache, il existe six promoteurs (1.1 ; 1.2 ; 1.3 ; 1.4 ; 1.5 et 2) qui dirigent la transcription du gène CYP19A1 dans les cellules de la granulosa. Le principal promoteur qui dirige la transcription au niveau de l’ovaire (cellules de granulosa) est le promoteur 2 (P2). Cependant, l’effet de la FSH et de l’IGF-1 sur l’activation de ces promoteurs d’aromatase demeure mal connu. De plus, la demi-vie du transcrit CYP19A1 est très courte avec une région 3’UTR relativement longue. L’analyse de la séquence 3’UTR montre la présence des motifs ARE (séquence riche en AU), des études antérieur montrent que ces séquences impliquent dans la régulation de la stabilité ou la dégradation de l’ARNm, ce qui est fort probable que la courte demi-vie de l’ARNm CYP19A1 est sous le contrôle post-transcriptionel. L’objectif de la thèse visait à étudier la régulation de l’expression du gène CYP19A1 chez la vache. Il y a deux thèmes soit étude de la régulation transcriptionnelle ciblant le promoteur et soit étude de la régulation post-transcriptionnelle impliquant la région 3’non traduite (3’UTR). Le premier objectif vise à étudier la régulation transcriptionnelle du gène CYP19A1. Nous avons étudié l'activité du promoteur ovarien bovin dans deux modèles de cellules de la granulosa, les cellules lutéinisées et nonlutéinisées in vitro, suite à une stimulation des cellules par la FSH ou IGF-1. Nous avons également évalué la voie de signalisation impliquée dans la régulation des différents promoteurs en utilisant un RT-PCR et un gène rapporteur (les différents promoteurs d’aromatase ont été insérés dans le vecteur pGL3promoter en amont du gène exprimant la luciférase). Les résultats de RT-PCR démontrent que la FSH et l’IGF-1 augmentent les concentrations d’ARNm provenant des deux promoteurs 2 et 1.1 dans les cellules de la granulosa non lutéinisées. Des expériences subséquentes ont montré que la FSH stimule le promoteur 2 via la voie PKA tandis que l'IGF-1 stimule le promoteur 2 via la voie PKC. La FSH et l’IGF-1 stimulent l’expression du promoteur 1.1 via la voie PI3K. L’analyse de l’activité luciférase démontre que dans les cellules de granulosa lutéinisées, la FSH stimule le promoteur 1.1 de façon dose dépendante et ne semble y avoir aucun effet significatif sur le promoteur 2. Nous avons donc comparé l’activité du promoteur PII/P2 humain, du rat, de la chèvre et de la vache dans les cellules de granulosa bovine lutéinisées. Le résultat le plus significatif est que le promoteur 2 bovine (et caprine) dépend de plusieurs facteurs de transcription (NR5A2, FOXL2) comparé au promoteur PII humain et celui du promoteur proximal du rat qui dépendent principalement de l'AMPc. En effet, nos résultats ont démontré une expression raisonnablement robuste du P2 bovine lorsque les cellules sont traitées à la forskoline, NR5A2 et FOXL2. Le facteur FOXL2 semble déterminer l'activité du promoteur 2 chez le ruminant. Le deuxième objectif vise à étudier la régulation post-transcriptionnelle du gène CYP19A1. Pour ce faire, nous avons déterminé la séquence minimale de l'ARNm CYP19A1 requise pour la régulation de sa demi-vie. Différents séquences de la région 3’UTR ont été insérés dans le vecteur pGL3promoter en aval du gène exprimant la luciférase ou soit dans le vecteur pGEMTeasy. Le vecteur pGL3promoter a été transfecté dans les cellules de granulosa lutéinisées pour évaluer l'impact de la séquence 3'UTR sur l'expression du gène rapporteur de la luciférase, alors que le vecteur pGEMTeasy a été utilisé pour la transcription in vitro afin de générer de l’ARNm. Ce dernier sera utilisé en réaction croisée au UV avec des extraits protéiques pour démontrer l’association du complexe ARNm/protéine. L’analyse de l’activité luciférase a permis d’identifier une séquence de 200 pb située entre 926 et 1134 pb de la région 3'UTR de l’ARNm CYP19A1 qui a réduit significativement l’activité de la luciférase. Selon les analyses de la réaction croisée au UV, une ou plusieurs protéines de 66 et 80 kDA se lient spécifiquement à la séquence de 200 pb qui réduit l’activité de luciférase. Cette protéine s'exprime dans les cellules de granulosa, mais n’a pas été détectée dans d'autres tissus comme le foie et le cœur. Par ailleurs, l’utilisation du gène rapporteur sensible à la FSH a suscité l’intérêt d'une compagnie pharmaceutique qui vend de l’equine chorionic gonadotropin (eCG) pour lui permettre de distinguer facilement l’eCG ayant une forte activité FSH et donc, avoir un produit commercial plus efficace et de meilleure qualité. Dans cette étude, nous avons développé un système de bioessai à la FSH basé sur la transfection des cellules avec un récepteur à la FSH et un gène rapporteur colorimètrique qui permet d’estimer l’activité de la FSH dans le sérum de la jument et qui pourrait être applicable au niveau de la ferme/industrie.

Genomic organization of human arylamine N-acetyltransferase Type I reveals alternative promoters that generate different 5'-UTR splice variants with altered translational activities

In humans, a polymorphic gene encodes the drug-metabolizing enzyme NATI (arylamine N-acetyltransferase Type 1), which is widely expressed throughout the body. While the protein-coding region of NATI is contained within a single exon, examination of the human EST (expressed sequence tag) database at the NCBI revealed the presence of nine separate exons, eight of which were located in the 5'non-coding region of NATI. Differential splicing produced at least eight unique mRNA isoforms that could be grouped according to the location of the first exon, which suggested that NATI expression occurs from three alternative promoters. Using RT (reverse transcriptase)-PCR, we identified one major transcript in various epithelial cells derived from different tissues. In contrast, multiple transcripts were observed in blood-derived cell lines (CEM, THP-1 and Jurkat), with a novel variant, not identified in the EST database, found in CEM cells only. The major splice variant increased gene expression 9-11-fold in a luciferase reporter assay, while the other isoforrns were similar or slightly greater than the control. We examined the upstream region of the most active splice variant in a promoter-reporter assay, and isolated a 257 bp sequence that produced maximal promoter activity. This sequence lacked a TATA box, but contained a consensus Sp1 site and a CAAT box, as well as several other putative transcription-factor-binding sites. Cell-specific expression of the different NATI transcripts may contribute to the variation in NATI activity in vivo.

The isolation, characterization, and application of WRKY genes as useful molecular markers in tropical trees

The improvement of tropical tree crops using conventional breeding methods faces challenges due to the length of time involved. Thus, like most crops, there is an effort to utilize molecular genetic markers in breeding programs to select for desirable agronomic traits. Known as marker assisted breeding or marker assisted selection, genetic markers associated with a phenotype of interest are used to screen and select material reducing the time necessary to evaluate candidates. As the focus of this research was improving disease resistance in tropical trees, the usefulness of the WRKY gene superfamily was investigated as candidates for generating useful molecular genetic markers. WRKY genes encode plant-specific transcriptional factors associated with regulating plants' responses to both biotic and abiotic stress. ^ One pair of degenerate primers amplified 48 WRKY gene fragments from three taxonomically distinct, economically important, tropical tree crop species: 18 from Theobroma cacao L., 21 from Cocos nucifera L. and 9 from Persea americana Mill. Several loci from each species were polymorphic because of single nucleotide substitutions present within a putative non-coding region of the loci. Capillary array electrophoresis-single strand conformational polymorphism (CAE-SSCP) mapped four WRKY loci onto a genetic linkage map of a T. cacao F2 population segregating for resistance to witches' broom disease. Additionally, PCR primers specific for four T. cacao loci successfully amplified WRKY loci from 15 members of the Byttneriae tribe. A method was devised to allow the reliable discrimination of alleles by CAE-SSCP using only the mobility assigned to the sample peaks. Once this method was validated, the diversity of three WRKY loci was evaluated in a germplasm collection of T. cacao . One locus displayed high diversity in the collection, with at least 18 alleles detected from mobility differences of the product peaks. The number of WRKY loci available within the genome, ease of isolation by degenerate PCR, codominant segregation demonstrated in the F2 population, and usefulness for screening germplasm collections and closely related wild species demonstrates that the WRKY superfamily of genes are excellent candidates for developing a number of genetic molecular markers for breeding purposes in tropical trees. ^

Microsatellite loci and the complete mitochondrial DNA sequence characterized through next generation sequencing and de novo genome assembly for the critically endangered orange-bellied parrot, Neophema chrysogaster

A suite of polymorphic microsatellite markers and the complete mitochondrial genome sequence was developed by next generation sequencing (NGS) for the critically endangered orange-bellied parrot, Neophema chrysogaster. A total of 14 polymorphic loci were identified and characterized using DNA extractions representing 40 individuals from Melaleuca, Tasmania, sampled in 2002. We observed moderate genetic variation across most loci (mean number of alleles per locus = 2.79; mean expected heterozygosity = 0.53) with no evidence of individual loci deviating significantly from Hardy-Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, and analyses indicated no evidence of null alleles across loci. De novo and reference-based genome assemblies performed using MIRA were used to assemble the N. chrysogaster mitochondrial genome sequence with mean coverage of 116-fold (range 89 to 142-fold). The mitochondrial genome consists of 18,034 base pairs, and a typical metazoan mitochondrial gene content consisting of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a single large non-coding region (control region). The arrangement of mitochondrial genes is also typical of Avian taxa. The annotation of the mitochondrial genome and the characterization of 14 microsatellite markers provide a valuable resource for future genetic monitoring of wild and captive N. chrysogaster populations. As found previously, NGS provides a rapid, low cost and reliable method for polymorphic nuclear genetic marker development and determining complete mitochondrial genome sequences when only a fraction of a genome is sequenced.

Structural studies of Helicase NS3 variants from Hepatitis C virus genotype 3 in virological sustained responder and non-responder patients

Background. About 130 million people are infected with the hepatitis C virus (HCV) worldwide, but effective treatment options are not yet available. One of the most promising targets for antiviral therapy is nonstructural protein 3 (NS3). To identify possible changes in the structure of NS3 associated with virological sustained response or non-response of patients, a model was constructed for each helicase NS3 protein coding sequence. From this, the goal was to verify the interaction between helicases variants and their ligands. Findings. Evidence was found that the NS3 helicase portion of non-responder patients contained substitutions in its ATP and RNA binding sites. K210E substitution can cause an imbalance in the distribution of loads, leading to a decrease in the number of ligations between the essential amino acids required for the hydrolysis of ATP. W501R substitution causes an imbalance in the distribution of loads, leading and forcing the RNA to interact with the amino acid Thr269, but not preventing binding of ribavirin inhibitor. Conclusions. Useful information is provided on the genetic profiling of the HCV genotype 3, specifically the coding region of the NS3 protein, improving our understanding of the viral genome and the regions of its protein catalytic site. © 2010 Rahal et al; licensee BioMed Central Ltd.